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Abstract 

We consider the problem of how to manipulate the entanglement properties 
of a general two-particle pure state, shared between Alice and Bob, by using 
only local operations at each end and classical communication between Alice 
and Bob. A method is developed in which this type of problem is found to 
be equivalent to a problem involving the cutting and pasting of certain shapes 
along with certain colouring problems. We consider two problems. Firstly, 
we find the most general way of manipulating the state to obtain maximally 
entangled states. After such a manipulation the entangled states |11) + 1 22} + 
. . . + | mm) are obtained with probability p m . We obtain an expression for the 
optimal average entanglement obtainable. Also, some results of Lo and Popescu 
pertaining to this problem are given simple geometric proofs. Secondly, we 
consider how to manipulate one two-particle entangled state to another 
with certainty. We derive Nielsen's theorem (which states a necessary and 
sufficient condition for this to be possible) using the method of areas. 



1 Introduction 

Quantum entanglement has many applications including quantum teleportation [p]], 
quantum cryptography 0, and quantum communication ||. This has led people to 
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regard entanglement as a resource. However, entanglement can exist in different forms 
and so it is useful to know how it can be manipulated from one form to another. In 
this paper the problem of manipulating a general pure two-particle entangled state in 
order to obtain maximally entangled states will be considered. We will also consider 
the problem of how to manipulate one general pure two-particle entangled state to 
another. Alice and Bob are allowed to do whatever they want locally and they 
are allowed to communicate classically with one another. They are not allowed to 
exchange quantum states. This type of situation has already been much discussed 
in the literature. The problem of how to optimally manipulate a large number, N, 
of copies of a general pure two-particle entangled states into maximally entangled 
states by local means has been completely solved in the asymptotic limit iV — > oo 
||. However, the perhaps more basic problem of how to manipulate a single copy 
of a general pure two-particle state into maximally entangled states has not been so 
extensively discussed. The most significant work on this is by Lo and Popescu [H 
who prove certain bounds relating to this problem. However their proofs, while being 
extraordinarily ingenious, are rather difficult to follow. The method developed in this 
paper, which completely solves the problem, involves the cutting and pasting of areas 
along with a colouring problem. Once the basic methods have been put in place, it 
is very easy to picture what is happening. This method is used to find the maximum 
obtainable average entanglement and to derive a formula of Lo and Popescu which 
gives the maximum probability of obtaining a given maximally entangled state. 

A related problem is how to transform one pure two-particle entangled state to 
another, and to establish which states will transform to one another in this way. 
Nielsen || has completely solved this problem. However his proof of his theorem uses 
some unfamiliar mathematics. An alternative reasonably simple proof of Nielsen's 
theorem is given here which, again, involves cutting and pasting of areas along with 
a certain colouring problem. 

2 Obtaining maximally entangled states 
2.1 Introduction 

The most general pure two-particle state can be written in the Schmidt form 
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where we choose Aj > Aj + i and where the states \i)a,b are orthonormal. We want to 
manipulate this state in order to obtain states which are of the form 

\<p m ) = —Y\k) A \k) B (2) 

We will call this state an m-state. An m-state is equivalent to log 2 m copies of 2- 
states Q. After the process is completed we should have a certain m-state with a 
certain probability p m . Particle A goes to Alice and particle B goes to Bob. Alice 
and Bob are allowed to perform whatever operations they want locally and also they 
communicate classically with each other. This can happen in the following way. Alice 
performs a measurement and communicates the result to Bob who then performs a 
measurement which depends on the result of Alice's measurement, and then Bob 
communicates his result back to Alice and she makes another measurement, and so 
on back and forth. This is most general way in which Alice and Bob can manipulate 
their state without actually exchanging quantum states in the process. Whatever 
measure of entanglement we employ, the amount of entanglement should not increase 
during such a process. Lo and Popescu show that, for the very special case of a 
two-particle pure state, this process is equivalent to one in which Alice makes one 
measurement and communicates the result to Bob who then may perform a unitary 
evolution operation on his particles. The reason for this significant simplification is 
that, due to the Schmidt decomposition, any operation by Bob is equivalent, so far as 
the resulting form of the state is concerned, to some operation by Alice. Hence, Alice 
can simply do everything herself in one go and then communicate the final result to 
Bob. The most general measurement Alice can make is a POVM. This is equivalent 
to Alice introducing an ancilla, S, performing a general unitary evolution on particle 
A and the ancilla S, and then making a projective measurement on the ancilla. Let 
us imagine that the ancilla has a basis set of states \l)s- Since we allow completely 
general evolution of A + S we can assume, without loss of generality, that the final 
measurement projects onto subspaces spanned by the states \l)s- Furthermore, it is 
shown in the appendix that there is no advantage to be had by performing a non- 
maximal (i.e. degenerate) measurement and so we can assume that this measurement 
is maximal and projects onto the operators \l)s{l\- For each outcome, I, the two 
particles A and B should be projected into an m-state |<^ m; ) ? where 

i mi 

Wm l Y = ^Y,\ k )A\k) l B (3) 

The superscript I is included since we do not require that the Schmidt vectors satisfy 
\k) l B = \k) l B for I ytz I', After Alice has communicated the result, I, of the measurement 
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to Bob, Bob could rotate these vectors into the same standard form for all I (thus 
removing the need for the superscript at this stage), but this is not important. It is 
enough that Alice and Bob know what I is so they know what state they have. We do 
not require a superscript / on the |/c)a states since, as explained below, Alice can rotate 
her Schmidt vectors to standard form as part of the overall unitary transformation 
she performs. Just before Alice makes her measurement projecting onto \l)s the state 
of the system will be 

I ^target) =EV/^IOsW (4) 
I 

where the coefficients y/JTi can be taken to be real since any phases can be absorbed 
into appropriately redefined \l)s- Note that if, at this stage, the \i)a states had a 
superscript / then they could be rotated into standard form by applying a series 
of controlled unitary operations, Ui, to A where the control is the \l)s state. At 
this stage Alice has not done anything which is irreversible. Having completed her 
local manipulations Alice will perform a maximal projective measurement. It follows 
from the result of Lo and Popescu that manipulating the state into the form (|J) and 
then measuring onto the basis is equivalent to the most general procedure for 
manipulating the two particles by local means to m-states. We will use equation 
dD later when we come to show that the method developed in the next section is 
equivalent to the most general method. 



2.2 How to obtain maximally entangled states 

Now consider the initial state \ip) given in equation ([]]). We will introduce an ancilla, 
R, in the state \1}r. This ancilla has basis states \ti)r where n = 1, 2, . . . N. We will 
take N to be very large and will want to consider the case where N tends to infinity. 
We define the integers iVj = NXi where, for the moment, we are taking the Aj to be 
rational numbers so the integers iVj can be found with N finite. This constraint can 
be relaxed when N tends to infinity. The initial state of RAB is jl)^^). Alice now 
evolves R + A using the following transformations. 

I Ni 

\l)n\i) A ^ ^=Q»*)|^ (5) 

V iV i n=l 

The fact that these transformations evolve orthogonal states to orthogonal states 
ensures that they can be implemented by unitary evolution. Under these transforma- 
tions the state |1)#|^>) evolves to 

l^start) = 12 \n)p.\i)A\i)B (6) 

n=l 
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Figure 1: The start state can be plotted on a graph. Such a graph is called an area 
diagram. This area diagram has step structure in which the steps go down towards 
the right. 



which we will call the start state. We see that each of the N terms in this superposition 
has the same amplitude. Each of these terms will be represented by a rectangular 
element width 1 and height 1/N. The area of the element is 1/N and equal to 
the probability associated with the corresponding term in (^j). Each element can be 
labelled by (n,i) corresponding to the term \n) A\j) b (initially i = j but after 
Alice has performed operations on her particles this will not necessarily be the case 
for every term). The elements are then arranged on a graph where elements having 
the same n are placed in the same row and elements having the same i are placed in 
the same column. The resulting graph looks like a series of steps as shown in Fig. |l[ 
We will call this the area diagram. The total area under the steps is 1 corresponding 
to the total probability. On the vertical axis n/N is plotted. On the horizontal axis 
the % which appears in \%)a is plotted. For small n every % position will be filled. 
However, because of the form of the state (J^j, once n gets bigger than Ni there will 
no longer be any \I)a\I)b terms. This is the reason for the step at n = Nj. There 
will be further steps at each n = iVj. The height of the rightmost step is Nj/N = Xj. 
Subsequent steps will be at heights A, as shown in Fig. [I]. If a projective measurement 
were to be performed on R at this stage then, corresponding to each outcome n, the 
state of A + B would be projected onto an m- state where m is equal to the number of 
elements in the nth row on the area diagram. This will give a distribution of m-states 
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with probabilities which are equal to the large horizontal rectangular areas of width 
m and height X m — A m+ i formed by extending the horizontal parts of the step back 
to the vertical axis. However, the area in the diagram can be moved around by Alice 
in a way to be described below by performing local unitary operations. When this is 
followed by a projective measurement on R different distributions of m-states can be 
realised. 

The terms \n) aU) b and \n') R \i') a\?) b are bi-orthogonal iff R(n\n') R A(i\i') a — 
(orthogonal at Alice's end) and b(3\j')b — (orthogonal at Bob's end). If two 
terms are only orthogonal at either Alice's end or at Bob's end then they are mono- 
orthogonal. In the area diagram we will impose the constraint that all elements in a 
row, that is for a given n, correspond to terms which are bi-orthogonal. This ensures 
that when we perform a projective measurement onto R the resulting state will be 
an rn-state. 

We will colour all the elements which correspond to terms which are mono- 
orthogonal to one another a given colour. Elements corresponding to terms which are 
bi-orthogonal will be coloured with different colours. Thus, initially, all the elements 
in a given column are the same colour and every column is coloured a different colour 
to every other column. When area is moved around the constraint that terms corre- 
sponding to elements in a row be bi-orthogonal means that all elements in any given 
row must be different colours. 

The method of areas to be developed here involves moving area elements around 
in such a way that there is a net movement of area up and to the left. We will see that 
it is possible to have a net movement of area up the area diagram but not down the 
diagram. This means that the net movement of area across any horizontal line drawn 
on the diagram must be up. The basic unitary operation, U(n,i <-> n',i'), employed 
by Alice is defined by the transformation equations 

\n) R \i) A -> \n') R \i') A (7) 

\n) R \i')A -> \n) R \i)A (8) 

with no change for all other \n") R \i")A- We will call this the swap operation. The 
effect of this operation is to move elements around on the area diagram. If there are 
elements at both the (n,i) and (n',i r ) positions then they will have their positions 
swapped. If, there is only an element at one of the two positions then it will be 
moved to the other position while the original position will become vacant. These 
moves will not effect Bob's part of the state. If two terms are bi-(mono)-orthogonal 
before the swap operation is applied to one or both of them then they will be bi- 
(mono)-orthogonal afterwards. In other words, the swap operation does not change 
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the colour of the elements. Initially, as stated above, all elements in the same row are 
different colours and elements in the same column are the same colour. In moving 
elements of area around we impose only the constraint that, at the end of the process, 
all elements in a row are different colours. Although elements in any given column 
start off being the same colour, we do not demand that this is true at the end of the 
process. We can move a large number of elements at once. In the limit as iV — > oo 
the elements will become infinitesimal in height. Hence, in this limit, we can make 
horizontal cuts anywhere. We can make vertical cuts along the edges of the columns. 
The area can be cut up into smaller pieces and then pieces can be moved around 
and pasted into new positions. It is possible to move area around like this in any 
way we want by repeated applications of the swap operation. The empty space above 
the steps can be used as a clipboard for the temporary storage of pieces of area to 
facilitate the rearrangement of area if required. 

In Fig. |2] the original step structure is shown by a full line. A new step structure is 
shown by a dashed line. We will impose the constraint (to be justified later) that the 
new step structure consists of steps which, like the original steps, go across and down 
(but never up) towards the right. The area S = Y^S r which will be cut away from 
some parts of the steps is equal to the area T = J2T r which will be added to other 
parts of the steps. Each of these areas consist of R smaller disjoint parts labelled by 
r going up the diagram. Note, area T r lies between area S r an SV+i- Note also that 
the areas S r and T r can themselves be made up of disjoint parts. The constraint that 
there is no net movement of area downwards means that YZ-=i T r < X£=i &r f° r au 
r'. The original steps are of height Aj. Let the new steps be of height A^. Since the 
columns are of unit width, these lengths are numerically equal to the areas of the 
columns, and hence the constraint that area is moved only to the left is equivalent to 
the set of constraints 

E K < E a, (9) 

i=p i=p 

for all r = 1 to I with equality holding when r = 1. 

We could simply move the area S to the area T element by element by applying 
the swap operator. However, if we did this it is likely that elements corresponding to 
mono-orthogonal terms would end up in the same row. If we can redistribute the area 
so that it corresponds to the new step structure but without there being any elements 
of the same colour in the same row then we will have realised another distribution 
of m-states. This is because Alice could then perform a projective measurement on 
R and corresponding to each outcome, n, will be a m-state and these m-states will 
clearly have a different distribution. The probability of a given m state is equal to the 
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Figure 2: Area is transferred up the area diagram to go from the start steps shown by 
the bold line to the target steps shown by the dashed line. The total area transferred 
up is S = T 



Figure 3: This figure shows the start steps (bold line) and the target steps(dashed 
line) of the reco louring problem discussed in the text. 

area of the horizontal rectangle formed by projecting leftwards the top and bottom 
of the step at position m. This rectangle has width m and height \' m — A' m+1 . Hence, 
for the new area diagram, the probability of getting an m state is 

Pm = (X' m - X' m+ i)m (10) 

We will show firstly that this colouring problem can be solved and secondly that the 
process described here is general in the sense that any distribution of m-states which 
can be achieved by local means can be achieved by this method. Hence, equations 
(H) and (|10"D define the possible distributions of m-states that can be obtained. 

The solution to the colouring problem will be explained by reference to the example 
shown in Fig. [3]. This example has R = 1 (since all of T is above all of 5"). However, 
it will be clear that the method works for the general case. We start by taking the 
rightmost column in the area S. This is area A in Fig |](a). This area is then swapped 
into the rightmost column of the area T so that point a coincides with point a' as 
shown in Fig. f|(b). However, it may be too long to all fit in to the right most column 
of area T (as is indeed the case in our example), and hence it will displace some of 
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(c) 



(d) 
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Figure 4: This figure shows the recolouring procedure described in the text. After 
the recolouring, no colour is in more than one place across any given row. 
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the area in the column below this, marked B in Fig. |](a). Because of the 

nature of the swapping operation, this area B will be moved to the old position of 
A so that point b coincides with point V as shown in Fig. ^(b). Area A is now in 
its final position. Let us imagine that the column from which area A was taken was 
coloured red. This red colour will now be divided between what is left of the original 
column and the area A in its new position. Since the steps go up towards the left, it 
is impossible to have red at more than one position in any given row. The rightmost 
column of T is now filled so we start on the next rightmost column. We swap area 
B into position in the next rightmost column of area T. Again, this could be too 
long. In our example it is too long and projects into the column below into the area 
marked C . Hence C will be moved to where B was. This takes us to the situation 
shown in Fig. f|(c). Now we move area C into position in the next rightmost column 
of area T. This could again be too long and project into the column below, but in 
this example this is not the case. Rather, area C is too short leaving a gap. Hence 
no area is moved back to the rightmost column of S and this means we have finally 
dispensed with the net effect of moving the area from this rightmost column of S. 
Now we select the next rightmost column of S. In our example this area is marked 
D in Fig. f|(c). This area is now moved up following the same procedure. Thus, we 
move D into position as high as possible in the rightmost column of T which has not 
yet been filled. This places it below area C. It could be the case that C is too long 
and projects into the column below in which case some area would be swapped back 
and we would have to continue as before. However, in our particular example area 
D fits in below C and the recolouring is finally completed as shown in Fig. |](d). 
This method can be applied to any recolouring problem of this nature. The general 
method is that area is swapped from the rightmost non-empty column of S to as 
high a position as possible in the rightmost non-full column of T. We note that a 
given colour can end up in, at most, two different columns and that any colour moved 
from a column will always end up higher than its original position. This means that 
it is impossible for two elements of area which started in the same column (having 
the same colour) to end up in the same row. Hence, the colouring problem has been 
solved. 

2.3 Proof that this is most general method 

Having shown that it is possible to have a net movement of area up the area diagram in 
any general way which is consistent with maintaining the step structure, we will now 
show that this corresponds to the most general way of manipulating entanglement to 
produce m-states in the sense that any distribution of m-states which can be achieved 



11 



by local operations classical communication can be achieved by this method. The idea 
of this proof will be to show that the target state (just before Alice measures on to 
the ancilla) can be put into a certain form which is inconsistent with any movement 
of area downwards. 

We have already established that the most general final state just before Alice 
makes her measurement (we call this the target state) is the state given in (|J). If we 
write Ui = fii/mi then @ becomes 



mi 



l«W> = T,EV^i\l)s\k}A\k) l B (11) 



I k=l 



We can consider further unitary transformations by Alice on this state to put it into 
a form in which every term has the same amplitude. Let the dimension of the ancilla 
S be M and define Mi = Mv\ (again we will let M — > oo) and let Alice perform the 
following transformations on (HI): 



|i>S--firr £ l»>S ( 12 ) 

M l n&W t 



where W t is the set of Mi integers from (Et=o M r) + 1 to E l r =o M r (we put M = 0). 
Under this transformation (0j) becomes 

i rni 

l*Uet> = 7^jEEE \n)s\k) A \k) l B (13) 
VA2 i k=in&Wi 

Every term in this state has the same amplitude. In arriving at this state from the 
previous target state (^JTT|) we have done nothing that is irreversible. Furthermore, 
if we measure onto the \n)s basis we are just as likely to get a given m-state as with 
the previous target state. Hence, we can regard this as our new target state. Any 
method by which m-states can be obtained is equivalent to manipulating the state 
into the form flT3"|). 



If we project onto \n)s we will obtain an m ra -state where m n can be read off 
from (0). We can relabel the n's such that m n+ i < m n . Thus, we can impose the 
following: 

Constraint A If outcome n corresponds to a m n -state then, without loss of gener- 
ality, we can impose the constraint that m n+1 < m n . 



By examining ( |T3"D we can see that there is a second constraint that can be imposed 
on the form of the final target state. 
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Constraint B For the target state we can, without any loss of generality, impose 
the constraint that 

\B(e\s(n\K**)\' < ^ (14) 
where \8)b is any normalised state for system B since this is true of equation 

We can identify the ancilla S with the ancilla R introduced earlier. Hence, S = R 
and M = N. Now consider the state ([I]). Since Alice's operations do not effect Bob's 
system we see that we have the following constraint: 

Constraint C While the state is being manipulated by local unitary operations 
by Alice, we will always have 

\B(m\ 2 = ^ = ^ as) 

for all j. 

Since (|5|) is related to (|T|) by reversible operations, we can take (^), which corresponds 
to an area diagram, as our starting point. We will now see that area cannot be moved 
down in the area diagram. For the purposes of this proof consider a change in the 
way the area diagram is plotted such that the j in \j)b (rather than the % in \%)a) 
is plotted on the horizontal axis. This will simply have the effect of redistributing 
elements horizontally but not vertically since n/N is still plotted on the vertical axis. 
We label elements in this modified area diagram by {n,j}. Constraint C implies that 
the column corresponding to a given j on this modified area diagram will always 
have the same area. This is because Alice's actions cannot effect the total area 
(or probability) associated with column j. However, Alice can change the n value 
associated with area elements and hence she can move area in column j up and down. 
It is possible that she can bring about a net movement of area (or probability) down 
in this column. This would lead to the area being compressed into a smaller space 
than it would "naturally" fit. Any net movement of area downwards, whether on the 
{n, j} picture or the (n,i) picture would correspond to this happening in at least one 
j column. This is exactly what we want to rule out. We will now see that any such 
net movement downwards will violate constraint B (which we were free to impose on 
the target state). 

If, at some stage, the state has been manipulated to the state then the area 
of the {n, j} element will be A n j = \B(j\s( n \^) | 2 - Initially, for the start state l^start) 
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in (0), we have 

n=l JV 

for all n! < iVj. However, if there is a net movement of area down the area diagram 
with respect to the horizontal line n = n' then, since the total area in a given column 
is constant (by constraint C), this net movement of area downwards must happen in 
at least one column of the modified area diagram. Hence, for at least one value of j, 
we must have 

n'<Ni / 

E < Set > ^ (17) 

n=l JV 

However, since M = N, equation (|I~4T ) implies 

n'<Ni / 

E <7 Ct < 77 (18) 

n=l iV 

The contradiction between (|T7| ) and ( JT8| ) proves that a net movement of area down- 
wards is not possible on the modified area diagram and hence neither is it possible 
on the unmodified area diagram. Note that this proof goes through for any sort of 
operations by Alice and in particular it does not assume that the only operations 
Alice can make are the swap operations defined in (|7|). 

To complete the proof that the manipulations described earlier are equivalent 
to the most general way of manipulating the state to obtain m-states we note the 
following. 

(i) The target state ([13]) can be represented on an area diagram in which n/N (where 

N = M) is on the vertical axis and m n is plotted on the horizontal axis. Since 
m n is an integer and since we can impose constraint A without loss of generality, 
this area diagram will have a step structure (in which the steps of integer width 
go down towards the right). 

(ii) The initial state can be taken to be the start state |^/ s t ar t) given in (j^) and this 

can be represented by an area diagram with the step structure. 

(iii) The total area of both these diagrams is 1. Therefore the most general way 
of manipulating the start state into m-states corresponds to going from one 
area diagram with the step structure to another with the step structure in 
a way consistent with the constraint that there is no net movement of area 
downwards. 
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(iv) The method, employing the swap operator, discussed previously can be used to 
go from one step structure to another in any way that is consistent with there 
being no net downward movement of area. Hence, it is equivalent to the most 
general method. 



2.4 Getting the highest possible average 

If we are only interested in the average amount of entanglement in the form of maxi- 
mally entangled states we can obtain, this being equal to E — J2 m Pm log 2 m, then it 
turns out that any movement of area will decrease this average. Hence this average 
has a maximum given by the original area diagram 

£ max = £(A m - A m+1 )mlog 2 m (19) 

m 

To see that any movement of area will decrease this average, consider moving one 
element on the area diagram (with area equal to -^) from the end of row A which has 
original width rriA, to the end of row B which has original width m B . Since we can 
only move area to the left we require that m B + 1 < m A (if m B + 1 = m A then the 
rows will simply have interchanged their lengths and hence there will be no change in 
the distribution of m-states). The original contribution of these two rows to E will 
be 

A.Ei n iti a i = — log 2 m A + -j-j- log 2 m B (20) 
The contribution afterwards will be 

A£ final = - 1) log 2 (m A - 1) + + 1) log 2 (m s + 1) (21) 
The difference between these two contributions is 



Affinal - AE initial = — log 2 



(m B + l\ ms fm A - l\ mA /m B + 1\ 
V m B J \ m A J \m A — lJ 



(22) 



It can only be advantageous to move elements of area if this quantity is positive. 
However, by making the substitutions 

1 1 

m A = x + - + r m B =x — - — r (23) 

we can see that A£fi na i— A£i mtia i is negative if r > 0. The constraint that m A > m B +l 
implies that r > and hence any movement of area must lead to a smaller E. We 
also see from this that since only one distribution of m-states leads to the maximum 
E, any attempt to alter the distribution of m-states will result in a decrease of E 
(and so we obtain another main result of Lo and Popescu). 
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Figure 5: This figure shows a strategy for obtaining the maximum probability of a 
given m-state. The area T is equal to the area S. This defines ro and A 
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2.5 Proof of a formula of Lo and Popescu 

We will now use this method to derive a formula central to the paper of Lo and 
Popescu [[|. Imagine that we have a general two-particle pure entangled state and 
we want to have a given m-state with as high a probability as possible. We want to 
know what this probability is and what strategy to use. This corresponds to the area 
redistribution shown in Fig. |5|. The target area diagram consists of a block of width 
m with an additional bit on top. This defines r and (see diagram). The area 

S has been moved to the area T where these two areas are equal. The height of the 
main block is A™!:* which can be calculated since we know that 

tit I o 



r A 



max 
m—ro 



U ro +T 



U ro + S 



E a, 



(24) 



i=m— ro+1 



where the last equality follows from the fact that the ith column is of area Aj. The 
total probability of getting the m-state is equal to the area of the main block, i.e. 



V 



" LA m-ro- 



We can use this formula if we know what r is since then A 



max 
m—ro 
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can be calculated from (^). This can be established by the following considerations. 
Define the area U r to consist of all the columns i = m — r + ltomin the start area 
diagram so that 

m 

*7r = £ ^ (25) 

i=m— r+1 

where r = 1, 2 . . . m. The area U r + T consists of a main block of width r and height 
A™^ x ro plus, for r 7^ r , an extra bit lying outside this block (when r < r a bit of T 
lies outside this block and when r > r a bit of 6+ lies outside the block). Hence, 

U r + T> rA™ x ro (26) 

with equality in the case r = Tq. Therefore, 

^ = min[i(C/ r + r)] (27) 

and we obtain the formula of Lo and Popescu 

171 

pST = nun-( E A,) (28) 

i=m— r+1 

where r = 1, 2, . . . m. The geometric origin of this formula is now clear. 



3 Proof of Nielsen's Theorem 
3.1 Introduction 

The set of constraints @ are exactly Nielsen's condition || for being able to ma- 
nipulate an entangled state, with Schmidt coefficients y/Xi to another, l^'), with 

Schmidt coefficients \fXi- However, we cannot immediately interpret the new area 
diagram as being equivalent to a state since, unlike in the original area diagram, 
a given column can be multicoloured. We will see that, nevertheless, we can use the 
area diagrams to prove Nielsen's theorem. This proof works along similar lines to the 
previous proof. First, we put the target state into step form. Alice introduces an 
ancilla S of dimension M with basis states \l)s- If the problem can be solved then, 
for similar reasons to before, she must be able to manipulate the total state into the 
form 

^targctHEv^lOsW (29) 
I 
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where 

WY -h^AWs (so) 

1=1 

Substituting ( p9|) into (|3~0"D we obtain 

l^targct) = E E V&|0sN)^N)B (31) 
i i 

Define Mu = /ijAjM. We apply the transformation 

\l)s\i)i - -y^r E \")s\i)A (32) 

where Wu is the set of integers from (Eto^ti) + 1 to J2k=o^ki (we set M 0i = 0). 
Under this transformation ( |31| ) becomes 

l*Uet> = 4jEE E |n>5|<>A|*>^ (33) 

V* 1 i neWu 

Each term has equal amplitude in this state. If we project it onto \n)s then we will 
get a state with, say, m n terms (where m n can be read off from (|33|)) not all of which 
are bi-orthogonal. We can relabel the n's so that m n > m n+ i and hence we can draw 
an area diagram with a step structure. The terms corresponding to a given n are not 
necessarily bi-orthogonal and hence we cannot impose on the target area diagram the 
constraint that elements in a row must all be coloured different colours. Rather we 
will have a different colouring problem. As before, we will identify the systems S and 
R so that M = N. 



3.2 Obtaining Nielsen's bound 

The start state is sta rt) in (0) and is represented by an area diagram with step 
structure in which each column is coloured with only one colour, this being different 
to the colour of the other columns. The target area diagram is shown in Fig. |6|. Since 
this has been recoloured the columns in this diagram will not necessarily be of one 
colour. The height of the ith column is = N[/N. We chose a number Q which 
we will let tend to infinity, though in the Fig. | we have set Q = 5. We divide each 
column up equally into Q pieces which are numbered q = 1, 2, . . . Q starting at the 
bottom. The gth piece in the ith column is labelled [q, i] If the target area diagram has 
been obtained from the start area diagram by moving finite sized bits of area around 
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Figure 6: The columns of the target diagram are divided up into Q equal parts. In 
this figure we have Q = 5 but we will let Q — > oo 
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(as will be the case in our recolouring strategy) then there will be a finite number of 
horizontal boundaries between different colours. Some of the pieces [q, i] are likely to 
have these boundaries on them. However, as Q — > oo the total area of such pieces will 
tend to zero and we can assume that each piece has a unique colour. The idea will be 
to collect all the pieces with a given q. Since their areas are proportional to they 
can correspond to the new state \ip') q . However, each of these pieces having the same 
q must be coloured with a different colour since the terms in \ip') q are bi-orthogonal. 
There are N- terms in the state vector corresponding to each column, and since these 
are divided into Q pieces, there are N-/Q terms corresponding to the [q, i] piece which 
will be of the form (unnormalised) 

1 qN'JQ 

^ q,i ^ = ~m ^ \n) R \i)A\j q i)B (34) 

V^ V n=(q-l)N'JQ+l 

The total state is the sum of all such terms. We can transform the total state by 
applying the transformation 

rQ qN'JQ 

\ j7j £ \n)n\i)A -»• \q)a\i)A (35) 

V n ={q-l)N'JQ+l 

for all q, i. This sends the term \[q,i\) to the state 

Fx- 

I [?,»]') = y -^\q)n\i)A\j q i)B (36) 

This transformation has simplified the state of the ancilla R for the terms correspond- 
ing to each piece [q, i\. In so doing we have recovered the coefficient \J-q- The total 
state is now 



W)-^HPM^a\3^b (37) 

V V 0=1 



Now, if we measure onto the basis we get a state which is a realisation of 
iff the terms are bi-orthogonal. In colouring terms, this means we require that all 
the pieces labelled with the same q in Fig. ^] should be of different colour. Thus, we 
have another colouring problem. If we can solve this colouring problem under the 
assumption that net movement of area is up, then we will have given a constructive 
proof that Nielsen's bound can be obtained. To complete the proof of Nielsen's 
theorem we need to show that area can only be moved up. This will be done later. 
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A way to solve this colouring problem was suggested to the author by A. Mahtani. 
This solution is obtained simply by correcting the solution to the previous colouring 
problem in Sec. 2.2 Firstly, this colouring procedure, illustrated by example in 
Fig. £|, is used to go from the start area diagram (representing \ip)) to the target 
diagram where the heights of the columns are A^. Now we note that it is a property 
of this reco louring procedure that a given colour can end up in two columns at most. 
There are two ways in which a piece of area can be moved to the left. Either it can 
be moved directly (as are the pieces A and D in Fig. [|, or it can first be swapped 
to the right and then be moved back to a further left position than it started in (as 
are pieces B and C in Fig. We will call the first type directly swapped pieces and 
the second type of pieces the swapped back pieces. Since these swapped back pieces 
(for example, piece B) must be shorter than the piece that displaced them from their 
original column (in the case of B this was piece A), since they always end up at the top 
of their destination column, and since the column they end up in is higher than their 
starting column, they must occupy proportionally less of their final column than the 
pieces that displaced them (piece B occupies proportionally less of its final column 
than piece A of its final column). This means that when the columns are divided 
up into a large number, Q, of equal pieces the colours of the swapped back pieces 
will not appear in two pieces with the same q. Hence, we will first correct the other 
colours. However, the procedure which corrects the colours of the directly swapped 
pieces disturbs this property of the swapped back pieces. Hence, after correcting for a 
bunch of directly swapped pieces, we will have to correct for the swapped back pieces 
as well. First we will consider the colours corresponding to the directly swapped 
pieces. Before we do that note that some columns will remain unchanged (neither 
have area moved into them, or out of them) and hence their colour cannot end up 
in two pieces with the same q. These columns can be completely ignored and we 
can consider only the remaining columns. We start at the rightmost column and go 
left considering only the columns which have changed. The first changed columns we 
meet will be shorter than originally, will each be coloured with only one colour, and 
will have had some area swapped out of them directly. After one or more of these 
monochromatic columns, we will meet a column into which these directly swapped 
pieces have been moved. The last directly swapped piece may displace a piece from 
this column which will be a swapped back piece S. This piece, when it is swapped 
back, will end up at a column, L, somewhere to the left. But before we get there, we 
may meet a few more monochromatic columns, Rk, from which pieces Xk have been 
directly swapped. We will count k — K, K — 1, . . . 1 backwards so first we meet Rk as 
we go leftwards. Eventually, on our journey leftwards, we meet column L into which 
S has been swapped, then pieces X K to X\. When the last piece Xi is swapped into 
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position it may displace piece S' which will be swapped back. Thus, we will repeat 
the same story as we continue leftwards. Columns will come in bunches of a few 
monochromatic columns (such as Rk) followed by a multicoloured column (like L). 
The columns L and Rk are shown in Fig. |7|(a) for the case K = 3. Columns not 
relevant to the present discussion are not shown. We label the original colour of L 
as Co and the colour of Rk as C&. Various distances (which are numerically equal to 
areas since the columns are of unit width) are marked on the figure. The strategy 
we will adopt is the following. We note that, as things stand, the colour Ck in all of 
Rk also appears in L and hence, there must be some of the same colour in different 
columns for the same q. To correct this, we can swap an area, of colour Co from L 
into Rk thus swapping the same area, of Ck back into L. We do this for all k. We 
then re-sort column L so that Yf. lies immediately below X^, both these areas being 
of the same colour Ck and so that colour C\ is above colour C2, etc. This is shown in 
Fig. |7|(b) (the role of area W will be explained later). We choose Yk to be such that 
the proportion of Ck in L relative to the height of L is equal to the proportion of Co 
in Rk relative to the height of Rk- This condition can be expressed as 

Xu -\- Yk Yu , . 

tk = k (3g) 

Furthermore, the area Yk of colour Cq is placed in Rk at the same relative position 
of Rk as the area X k + Y k of colour Ck has been placed in column L (see Fig. 0(b)). 
This ensures that when the columns are divided up as shown in Fig. |^, but with Q 
large, the pieces [q, i{\ in L of colour Ck will have a different colour from the pieces 
[q, iR k ] in Rk since the latter pieces will be of colour Co. 

Having carried out this correcting procedure for these columns we see that there 
is a problem. Piece S' which has been swapped back somewhere to the left of L is 
of colour Co- This piece may overlap (in the sense of occupying some pieces with the 
same q) with the piece Yi, also of colour Co, which is in column Ri. We can see S' will 
not also overlap with the pieces Y 2 , Y$, . . . (also of colour Co) in columns R2, R3, ■ ■ ■ 
for the following reasons: (1) It is smaller than X± and hence can only partly overlap 
with the piece X\ + Y\ (of colour Ci) in column L in Fig. ^ (2) This piece of colour 
C\ in column L does not overlap with the colour Co in columns R2, R3, ■ ■ ■■ Hence, 
this problem only concerns piece S' and column Ri . Let us assume that piece S' is in 
column V and that the original colour of column V is C . We can use colour C' Q to 
correct for piece S' in V and Y\ in R\ by essentially the same correcting procedure as 
before. Thus, we swap a piece of C and of area W from L' to into Y\ in R\ and at 
the same time swap a piece of colour Co and of area W from Y\ in Ri into V . The 
areas 5" and W are collected together at the top of column L. The size of area W 
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Figure 7: To solve the colouring problem we consider correct bunches of columns 
such as those shown here. We correct the diagram by moving the areas 



is chosen to be just such that S' + W (which is of colour C ) no longer overlaps with 
any of the colour Cq in Ri. The maximum size of area W is given by 

g + ^max _ ^max 

We could have W smaller than the value given by this equation since it is possible 
that not all of S' overlaps with Yi in Ri. Hence, 

W < S \ (40) 

This correcting procedure is carried out in the following way. First the rightmost 
bunch of columns like Rk and L are selected. Then the directly swapped pieces are 
corrected. And then the swapped back pieces are corrected. Then the next bunch of 
columns are subject to the same correcting protocol until all the bunches have been 
corrected. If this procedure can be carried out successfully then we will have solved 
this colouring problem. The only possible problem would be if we had to swap more 
of the colour Cq from L than there is in the column. However, we can show that this 
will not happen. From the Fig. |7|(a) we see that the area, z, of the original colour Co 
in L satisfies 

X' L -Y,X k -S = z>X' Ri >X' Rk (41) 
k=\ 

The inequality must be satisfied since otherwise elements of colour C\ in row L will 
be in the same row as elements of colour C\ in row R\ but we proved in Sec. II. 
B. that this could not happen with this colouring of the diagram. From ( |38| , f4TD we 
obtain 

Y k = - \ X k < X ' R * Xk < ^ (42) 

X' L — X Rk J2k> Xk> + S J2k> Xk' + S 

Since 5" had to be corrected for, S must also be corrected for (if L and Rk are taken 
to represent a general bunch of columns). To correct for S will require an area W of 
colour Co to be swapped into a column R'{ somewhere on the right. By analogy with 
( |40| ) we have 

W < - — 5l- < ^± (43) 

\' L — \ R n J2k> Xk' + S 

where the second inequality follows since \' R „ < < z. Hence, 

Y, Y k + W<z (44) 
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Figure 8: The mini-step form consist of the pieces being arranged along the i-axis 
starting with the q — 1 pieces. 



Which means that there is enough of the colour Co in column L to complete this 
colouring strategy. Hence, the colouring problem has been solved. 

3.3 Proof that Nielsen's bound cannot be beaten 

We need to prove that Nielsen's bound cannot be beaten. This is equivalent to proving 
that it is not possible to have net movement of area down the area diagram. The start 
state corresponds to a step structure. Each column in this can be divided up into 
Q pieces as shown in Fig. || Next we collect together all the pieces corresponding 
to a given q and place them in order along the i axis as shown in Fig. [§]. Thus, 
going along the i axis we have the q = 1 pieces, then the q = 2 pieces and so on. 
This can be accomplished by Alice performing swap operations. We will call this the 
mini-step form for the area diagram. (We are, of course, assuming that system A has 
a large enough Hilbert space to be able to do this. If this is not the case then an 
additional ancilla could be introduced to effectively increase the size of A's Hilbert 
space.) Let the state corresponding to this diagram be |^™art)- Now we change from 
the (n, i) to the {n, j} picture where j is the j in \j)b- The {n, j} element can only be 
moved up and down the diagram when Alice performs local operations. Hence, in the 
{n,j} picture all the mini-steps will overlay each other so there will be Q elements 
at each position. The height of the jth column will be Nj/Q. Hence, if we define 
A nj = | B (j| s (n|^ mini )| 2 , then for the start state we have 

n ' n'Q 

n=l iV 
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Figure 9: This shows an intermediate state of the area diagram involved in proving 
that Nielsen's bound cannot be beaten. 



for all n' < Nj/Q. The factor Q comes from the fact that there are Q sets of mini-steps 
overlaying each other in the {n,j} picture. 

Now we go back to the (n,i) picture and consider the target state fl29|). We can 
apply the transformation 

105 - 4t E \Q)s (46) 



V v l q£W l 

where Vi = jj[Q, Wi is the set of integers from (Xt=o K-) + 1 to Y/ r =o K ( we se t Vo = 0). 
We will let Q — ► oo. Under this transformation (EUI) becomes 



I target 



77 E \l)sW) q 



(47) 



where the superscript q is equal to I for q 6 Wj. This state now consists of a number 
of terms \q)sW) q each having the same amplitude Each term, \q)s\i J ') q , can 
individually be put into the step form 



\q)8\W^\V,q) 



gN/Q 

E 

i n=(q-l)N/Q+l 



E E \ n )s\i)A\i) Q 



B 



(48 



by applying transformations similar to If this transformation is applied to all 
terms then the resulting area diagram will consist of a series of mini-steps lined up 
vertically as shown in Fig. |9|. Next, Alice applies swap operations to move these sets 
of mini-steps so that they are lined up along the i axis starting with the q = 1 pieces 
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giving an area diagram in mini-step form (like in Fig. The state becomes 

rr n/q 

l*££t> = \//vEEE \n)s\j + I(q-l))A\j) q B (49) 

' q j n=l 

For this state we have 

\B(e\ s (n\^ t )\ 2 < ^Y, 1 = ^ (so) 

for any normalised state \6) B . 

The problem is to go from the start diagram in mini-step form to the target 
diagram which is also in mini-step form. If and only if we can do this can we also go 
between the corresponding diagrams in standard step form since Alice can transform 
reversibly between the two types of form of the area diagram. If there is to be net 
movement of area downwards in the mini-step form then this must happen for at least 
one value of j. Hence, comparing with (^3j), net downward movement of area implies 

fx? et >^ (5i) 

n=l iV 

for at least one value of j and n' < Nj/Q. However, (|50| ) implies 

E <T = E \B(j\s(n\*Ztt)\ 2 < ^ (52) 

n=X n=l iV 

which contradicts (EH]) and hence there can be no net movement of area downwards 
in the mini-step form. The standard step form area diagrams are simply elongated 
versions of one set of mini-steps in the mini-step form, and hence, by the similarity of 
these shapes, there can be no movement of area downwards in the standard picture. 
This proves Nielsen's bound (given algebraically in @). 



4 Conclusions 

In this paper a method of areas has been developed which enables us to understand 
the manipulation of pure two-particle entanglement. This approach has been used 
to find the most general way of transforming a general two-particle pure state into 
maximally entangled states. Certain results of Lo and Popescu were given geometric 
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interpretations. This method has also been used to prove Nielsen's theorem which 
pertains to going from one two-particle pure state to another with certainty. There 
remain a number of open problems relating to manipulation of two-particle pure 
entanglement which it may be possible to solve using the method of areas. Firstly, 
we could generalise Nielsen's theorem to the problem where we go from one state to 
another but not necessarily with certainty. Secondly, we could consider the problem 
of going from one state to a distribution of states. The method may also generalise 
to more than two particles (though it is not presently clear how this generalisation 
will work). 
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Appendix 

In this appendix we show that there can be no advantage if Alice makes a non- 
maximal rather than a maximal measurement onto S. Assume that the state just 
before measurement is 

5>IOs|*i> ( 53 ) 
i 

where is some state of system AB and not necessarily an m-state. Imagine 
that the projective measurement is non-maximal and one of its projectors is 1 1)5(1 1 + 
|2)5(2|. In the case of having the corresponding outcome, the resulting (unnormalised) 
state will be 

c 1 |l)s|$ 1 ) + c 2 |2) s |<l> 2 ) (54) 

This could, for example, be an m-state if system S is regarded as being part of sys- 
tem A. Rather than performing this non-maximal measurement Alice could instead 
change her notation for the \i)a states such that, if \i)a appears in the expansions of 
|$i) and I $2) she writes \i)a as \1,i)a- The remaining vectors \i)a are relabelled as 
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1 2, i)a- We are free to assume that the dimension of A is big enough to do this. Then 
we write \k,i)^ = Now Alice performs the transformations 

\1)s\1)a>^\2)s\1)a> (55) 

\2) s \1)a> -> \2)s\2)a> (56) 
Under these transformations the first two terms in ( |53"D become 

|2) s (c 1 |l)^# 1 ) + c 2 |2) A ,|$ 2 )) (57) 

A maximal measurement will now give rise to a state with the same form as the 
state as in ([54"D for the outcome 2. This trick can be repeated everywhere there is a 
degeneracy in the original non-maximal measurement and a maximal measurement 
can then be performed instead. This maximal measurement will give rise to the same 
distribution of the same states as the non-maximal measurement and so there can be 
no advantage to performing non-maximal measurements. 
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